subroutine load_s(ie, fij, fab, fia, & 
          iijkl, iijka, iabij, iiabj, icore, escf) 
!use common
use mod_iop
use mod_size
use mod_orbit
use mod_ioff 
implicit none
real*8, intent(inout) :: escf 
real*4, intent(inout) :: ie(*)
real*4, intent(inout) :: fia(*), fab(*), fij(*), iijkl(*)
real*4, intent(inout) :: iijka(*), iabij(*)
real*4, intent(inout) :: iiabj(*), icore(*)
!
integer nsize, iend, iaabb, i, j
integer iaaaa, ibbbb, iabab, ibaba, iabba, ibaab 
integer, external :: idsymsoc, irpdso
real*8, external :: dnrm2
logical yesno
!
open(unit=77, file='energydataus', form='unformatted')
rewind 77
read(77)(ie(i), i = 1, noa+nva)
read(77)(ie(i), i = noa+nva+1, noa+nva+nob+nvb)
read(77) ESCF
close(77)
!
call zeroso(fab, nifvv)
call zeroso(fij, nifoo)
call zeroso(fia, nifvo)
inquire(file='FAIAAS__', exist=yesno)
if(yesno) then 
call getallso(fia(ifvo(1)), 1, 'FAIAAS__')
do i = 1, nva
!   write(6,*)'fai', i, fia(ifvo(1)-1+i)
enddo

call getallso(fia(ifvo(2)), 1, 'FAIBBS__')
if(ldebug) then 
call printdnrm2('fai', fia, 22)
endif
endif
!
!======IJKL======
!
iaaaa = 1
ibbbb = iaaaa + idsymsoc(1, popa, popa, 1, popa, popa, 1) 
iabab = ibbbb + idsymsoc(1, popb, popb, 1, popb, popb, 1)
iend  = iabab + idsymsoc(1, popa, popb, 1, popa, popb, 1)
!
call getallso_s(iijkl(ioi1(1)), 1, 'IJKLAAS_')
call getallso_s(iijkl(ioi1(2)), 1, 'IJKLABS_')
call getallso_s(iijkl(ioi1(3)), 1, 'IJKLBBS_') 
if(ldebug) then 
call printdnrm2('iijkl', iijkl, 1)
endif
!
!======IJKA======
!
call getallso_s(iijka(ioi2(1)), 1, 'IJKAAAS_')
call getallso_s(iijka(ioi2(2)), 1, 'IJKAABS_')
call getallso_s(iijka(ioi2(3)), 1, 'IJKABBS_')
call getallso_s(iijka(ioi2(4)), 1, 'IJKABAS_')
if(ldebug) then 
call printdnrm2('iijka', iijka, 2)
endif
!
!======IABIJ======
!
call getallso_s(iabij(ioi3(1)), 1, 'ABIJAAS_')
call getallso_s(iabij(ioi3(2)), 1, 'ABIJABS_')
call getallso_s(iabij(ioi3(3)), 1, 'ABIJBBS_')
if(ldebug) then 
call printdnrm2('iabij', iabij, 3)
endif
!
!======IIABJ======
! IIABJ<i,a|b,j> would be sort as (bi,aj)
!
if(iowmbej==0) then 
  call getallso_s(iiabj(ioi4(1)), 1, 'BIAJAAAS') 
  call getallso_s(iiabj(ioi4(2)), 1, 'BIAJABAS') 
  call getallso_s(iiabj(ioi4(3)), 1, 'BIAJBBBS')
  call getallso_s(iiabj(ioi4(4)), 1, 'BIAJBABS')
  call getallso_s(iiabj(ioi4(5)), 1, 'BIAJAABS')
  call itranspso(iiabj(ioi4(5)), iiabj(ioi4(6)), & 
     vrta, popa, 1, vrtb, popb, 1, 1)
elseif(iowmbej==1) then 
  call clmoioso0_s(1, icore, vrta, popa, 1, vrta, popa, 1, 'WMBEJAAA')
  call clmoioso0_s(1, icore, vrtb, popb, 1, vrtb, popb, 1, 'WMBEJBBB')
  call clmoioso0_s(1, icore, vrta, popb, 1, vrta, popb, 1, 'WMBEJABA')
  call clmoioso0_s(1, icore, vrtb, popa, 1, vrtb, popa, 1, 'WMBEJBAB')
  call clmoioso0_s(1, icore, vrta, popa, 1, vrtb, popb, 1, 'WMBEJAAB')
  call clmoioso0_s(1, icore, vrtb, popb, 1, vrta, popa, 1, 'WMBEJBBA')
else 
  write(6,*)'error of iowmbej'
  stop
endif
!
!======IABCI======
!
return
end 
